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THz  Imaging 

of  Carrier  Concentration  and  Mobility  in  Silicon. 

Michael  Herrmann,  Masahiko  Tani,  Kiyomi  Sakai,  Ryoichi  Fukasawa 


Abstract  - Silicon  samples  have  been  imaged  with  THz 
radiation.  The  samples  are  n-type  substrates,  with  half  of 
one  surface  ion  implanted  with  boron  resulting  in  a 0.55  pm 
thin  strongly  p-doped  layer.  Silicon  reflects  and  absorbs  the 
THz  beam  due  to  a plasma  resonance,  subsequently  the 
additional  layer  results  in  strong  image  contrast  depending 
on  the  implantation  dosage.  Spectra  of  these  samples  have 
been  fit  using  a refractive  index  and  extinction  coefficient 
calculated  from  the  plasmon  resonance:  carrier  concentrat- 
ion and  relaxation  time  are  used  as  adjustable  parameters. 
By  the  fitting  spectra  at  every  image  position,  THz  images  of 
the  carrier  concentration  and  carrier  mobility  in  silicon  can 
be  obtained. 


I.  Introduction 


II.  Experimental 

The  optical  setup  used  for  the  experiments  is  shown  in 
Fig.  1.  The  output  of  a Ti: Sapphire  laser  is  split  into  a 
pump  and  probe  beam.  The  pump  beam  hits  the  THz 
emitter,  a dipole  antenna  on  LT-GaAs.  The  emitted  THz 
radiation  is  focused  with  parabolic  mirrors  on  the  sample, 
which  is  scanned  in  two  directions  vertical  to  the  beam. 
The  transmitted  radiation  is  focused  on  the  THz  receiver, 
another  dipole  antenna,  which  is  gated  by  the  probe 
beam.  When  a THz  pulse  and  a probe  pulse  arrive  at  the 
receiver  simultaneously,  a signal  is  generated  and 
recorded  with  a lock-in  amplifier.  The  probe  pulse  can  be 
delayed  with  a delay  stage.  In  this  way  a time-domain 
waveform  is  recorded  at  every  image  position. 


Silicon  is  one  of  the  technologically  most  important 
materials  today.  It  is  also  one  of  the  most  cleanly 
produced  materials,  and  consequently  the  characterisation 
of  silicon  is  an  essential  part  of  the  silicon  production 
process  securing  the  quality  of  a delicate  product.  With 
the  plasmon  frequency  of  silicon  in  the  THz  regime,  THz 
radiation  should  be  suitable  for  optically  determining  the 
electrical  properties  of  silicon.  Some  recent  studies  on 
silicon  and  other  semiconductor  spectra  in  the  THz 
regime  can  be  found  in  Refs.  1-5. 

In  1995  it  was  demonstrated  that  THz  frequency  beams 
can  be  used  in  combination  with  a scanning  technique 
providing  local  resolution  of  samples  [6  - 8].  THz 
Imaging  (also  called  "T-ray  Imaging")  may  among  other 
prospective  applications  serve  as  a convenient  tool  for  the 
non-destructive,  contactless  and  locally  resolved  charac- 
terisation of  silicon  and  other  semiconductors.  The  use  of 
THz  Imaging  for  semiconductor  characterisation  was  first 
proposed  by  Mittleman  etal  [8,  9]  in  1996.  They  used  a 
specialised  THz  Imaging  setup  with  a magnetic  field  and 
a double  detector  system  that  makes  use  of  the  rotation  of 
the  polarisation  of  THz  light  due  to  the  "terahertz  Hall 
effect". 


The  samples  are  three  specimen  cut  from  a 380  pm  thick, 
phosphorus-doped  n-type  Si  wafer  with  4.5  f2cm 
resistivity.  Half  of  one  surface  of  each  sample  has  a 0.55 
pm  thin  p+-type  layer  produced  by  ion-implanting  boron 
at  a dosage  of  5xl013  cm*2,  5xl014  cm'2  or  5xl015  cm*2, 
respectively,  and  subsequent  annealing. 


In  this  paper,  we  propose  the  characterisation  of  silicon 
with  a standard  THz  Imaging  setup.  We  investigate  n- 
and  p-type  silicon,  and  a "thick"  wafer  as  well  as  thin  ion- 
implanted  layers.  The  measurements  are  compared  to 
analytic  and  computer-supported  calculations. 
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Nishi-ku,  Kobe  651-2401,  Japan 
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Fig.  2:  Schematic  of  the  sample.  Half  a surface  has  been  ion- 
implanted  with  boron  resulting  in  a p+-type  layer. 
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III.  Intensity-Based  Images 


In  the  intensity-based  THz  images,  the  non-implanted 
sides  of  the  samples  transmit  between  40  to  50  % of  the 
THz  signal  (i.e.  electric  field),  and  the  implanted  layers 
provide  another  drop  in  signal  intensity  that  depends  on 
the  doping  concentration.  Fig.  3.  shows  a THz  image  of 
the  5x10l  cm'2  Si  sample.  On  part  of  the  image  the  THz 
beam  was  allowed  to  pass  the  sample.  The  corresponding 
(right-hand)  part  of  the  image  is  bright,  representing 
maximum  transmittance,  or  the  reference  level.  The  rest 
of  the  image  shows  two  regions  of  different  transmittance 
corresponding  to  the  implanted  (top)  and  the  non- 
implanted  (bottom)  part  of  the  sample.  On  all  intensity- 
based  images  shown  here,  zero  brightness  (i.e.  black) 
represents  zero  transmittance. 


Fig.  3:  THz  image  of  the 
5xlOI:i  cm'2  Si  sample. 
On  the  right-hand  side 
the  THz  beam  reaches 
the  receiver  directly,  re- 
sulting in  a strong  signal. 
The  sample  (larger,  left 
part  of  the  image) 
reflects  and  absorbs  the 
THz  beam  partly.  The 
upper  part  of  the  sample 
is  the  implanted  side. 
Image  size  is  20  mm  x 
20  mm.  The  dark  areas  in 
the  edges  are  shadows 
from  the  sample  holder. 

In  Fig.  4 THz  images  of  the  three  samples  are  compared. 
Their  undoped  sides  look  identical,  but  the  doped  sides 
differ  according  to  doping  concentration.  The  weakly 
implanted  layer  causes  only  a slight  reduction  in 
intensity,  while  the  strongly  doped  layer  reduces  the  THz 
signal  by  another  factor  5.  This  is  shown  quantitatively  in 
Fig.  5 which  is  based  on  vertical  intersections  through  a 
set  of  images  as  those  in  Fig.  4. 

IV.  Theory 

The  interaction  between  a semiconductor  and  electro- 
magnetic radiation  at  THz  frequencies  is  due  to  plasmon 
resonance  and  can  be  described  by 

s((o)  = £oa P --  - n 2 = (n  + ik)  (1) 

col  -icor  1 

[10]  where  0)p  is  the  plasmon  frequency  with 


Fig.  4:  Intensity-based  THz  images  of  three  n-type  Si  samples 
the  lower  parts  of  which  have  been  ion-implanted  with  boron, 
resulting  in  0.55  pm  thin  p+-type  layers.  Implantation  dosages 
(from  left  to  right):  5xl013  cm'2,  5xl014  cm'2  and  5xl015  cm'2. 
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Fig.  5:  Intersections  through  THz  images  of  the  three  samples. 
While  the  92  % of  the  THz  electric  field  signal  can  transmit  the 
weakly  doped  layer  (dotted  line),  only  18  % pass  through  the 
strongly  doped  layer  (solid  line). 
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scattering  time,  e (or  £«>)  the  dielectric  function  (at  infinite 
frequency,  respectively),  n = n + ik  the  complex 
refractive  index,  and  N,  //,  m*  e and  e0  the  carrier 
concentration,  mobility,  effective  mass,  electron  charge 
and  vacuum  permittivity,  respectively. 

If  the  semiconductor  properties  are  known,  the  complex 
refractive  index  can  be  calculated  from  Eq.  1,  and 
transmittance  through  a single  semiconductor  layer  can 
then  be  calculated  using 


rp  . i{n-\)-d 

tU 4n  x e ‘ (4) 

^ref  (n  +1)  (n-l)Z  2m—  d 

(«+l)2 

[12]  where  the  middle  fraction  describes  reflection  at  the 
surfaces  and  in  the  right  fraction  the  numerator  represents 
bulk  absorption  and  phase  shift  and  the  denominator 
takes  account  of  multiple  reflection. 

For  the  implanted  sides,  a two-layer  model  has  to  be  used 
as  described  by  Eq.  5. 
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[12].  However,  much  of  what  happens  in  the  layer  can 
also  be  understood  from  a thin-layer  approximation  of 
Eq.  4.  For  a given,  constant  implantation  dosage  Nar,  the 
carrier  concentration  N = Na/d  increases  with  decreasing 
thickness  d , so  that  h2  is  proportional  to  1/d  in  Eqs.  1 and 
2.  Using  this  and  going  towards  zero  thickness,  Eq.  4 
works  out  to 
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In  a strongly  doped  sample  r can  drop  to  10  fs  and  is  then 
« l/co  at  THz  frequencies,  so  that  the  imaginary  part  of 
Eq.  6 can  be  neglected  in  determining  the  absolute  value 
of  the  transmittance.  The  phase  shift  can  be  calculated 
from  Eq.  6 to  be 


tan  <p(a> ) = 


ImfrrtwQp)) 

XsktUni**)) 


(7) 


in  a further  approximation  valid  for  Nar » 10 14  cm2  and 
thus  at  least  for  the  strongly  doped  sample.  Using  these 
approximations  the  delay  time  in  a thin  layer  equals 
minus  the  relaxation  time,  At  = (p(co)/co  = -r. 


V.  Fitting  Spectra 

Fig.  6 shows  a three-parameter  fit  of  THz  transmittance 
for  the  non-implanted  side  using  Eqs.  1 to  4.  With  1960 
cm2/Vs  the  resulting  mobility  is  significantly  higher  than 
the  literature  value  of  1450  cm2/Vs  for  intrinsic  silicon. 
They  are  in  fairly  good  agreement  with  those  of  van  Exter 
et  al.  [1,  2]  who  also  measure  significantly  increased 
values  for  (u . From  the  resulting  values  for  N and  //,  the 
resistivity  is  calculated  to  p = 4.4  Qcm  which  is  in  very 
good  agreement  with  a value  of  p = 4.5  Qcm  measured 
on  this  sample  with  spreading  resistance  technique. 


Fitting  n-typc  Si  substrate  with  onc-laycr  mode! 


Frequency  (THz) 


Fig.  6:  A THz  spectrum  measured  by  Fourier  Transform 
Infrared  Spectrometry  (FTIR)  of  the  non-implanted  side  was  fit 
using  Eqs.  1 and  2. 


transforming  the  waveforms  and  calculating  the  time 
delay  from  the  FFT  phase  shift  at  various  frequencies. 
Fig.  7 shows  a THz  image  based  on  time  delay  at  1.25 
THz  determined  in  this  way.  The  delay  of  the  implanted 
side  relative  to  the  non-implanted  side  amounts  to  -8.5  ± 
2 fs.  The  calculated  value  is  -6  fs.  The  calculation  used 
Eq.  4 without  the  correction  for  multiple  reflection  for  the 
non-implanted  side  and  Eq.  5 for  the  implanted  side. 
Multiple  reflection  in  the  substrate  was  excluded  from  the 
formulas,  since  the  recorded  time-domain  waveforms 
were  short  and  did  not  include  pulses  reflected  through 
the  substrate.  Multiple  reflection  within  the  doped  layer 
however  was  included. 


Fig.  7:  Time-delay  THz  image  (17  mm  x 5 mm)  of  the  strongly 
doped  sample.  The  implanted  region  at  the  left  appears  dark 
representing  earlier  arrival  time  by  8.5  ± 2 fs.  The  dark  vertical 
line  in  the  middle  is  probably  a diffraction  effect. 


VII.  Conclusion 

We  have  used  a conventional  THz  Imaging  system  for 
investigating  the  local  distribution  of  carrier  concen- 
tration and  mobility  in  silicon  samples,  n-  and  p-type 
silicon  have  been  measured,  and  a "thick"  wafer  as  well 
as  thin  ion-implanted  layers.  The  measurements  have 
been  compared  to  analytic  and  computer-supported 
calculations.  For  a thick  sample  (i.e.  > several  10  pm), 
both  carrier  concentration  and  mobility  can  be 
determined  by  fitting  transmission  spectra.  Mobility  data 
is  significantly  increased  against  the  literature  value,  but 
we  see  excellent  agreement  with  the  nominal  value  of 
sample  resistivity.  For  a thin  layer,  either  concentration  or 
mobility  can  be  determined  in  this  way,  the  other  quantity 
being  accessible  via  time-delay  information.  The  time- 
delay  of  a 0.55  pm  thin  p+-type  layer  on  n-type  substrate 
has  been  determined  to  -8.5  ± 2 fs  in  fair  agreement  with 
calculations  using  a two-layer  model.  The  method  is  non- 
destructive, does  not  need  electric  sample  contacts  and 
has  a local  lateral  resolution  of  about  1 mm. 


In  the  thin  layer  approximation,  the  transmittance 
depends  only  on  the  product  N/j,  not  on  N and  jj. 
separately  (cf.  Eqs.  3 and  6 with  the  imaginary  part  in  Eq. 
6 neglected).  This  is  why  only  one  of  these  two  quantities 
can  be  fit  from  a transmittance  spectrum.  In  order  to  fit 
the  second  quantity,  phase  information  must  be  taken  into 
account. 

VI.  Time  Delay 

Since  the  time  delay  is  much  smaller  than  the  pulse 
duration,  determining  the  time  delay  from  specific 
features  of  the  measured  THz  waveforms  is  questionable 
because  even  slight  changes  of  pulse  shape  due  to 
dispersion  can  change  the  result  significantly.  However, 
relatively  good  agreement  is  reached  by  Fourier 
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